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Unstable Semiclassical Trajectories in Tunneling 
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Some tunneling phenomena are described, in the semiclassical approximation, by unstable com- 
plex trajectories. We develop a systematic procedure to stabilize the trajectories and to calculate 
the tunneling probability, including both the suppression exponent and prefactor. We find that 
the instability of tunneling solutions modifies the power-law dependence of the prefactor on ft as 
compared to the case of stable solutions. 



Tunneling in systems with many degrees of freedom 
has been a subject of continuous theoretical research for 
the last decades The interest is heated up by the 
recent experimental observations 0, H, 01 of non-trivial 
dynamical properties of multidimensional tunneling. In 
the semiclassical framework tunneling is described as mo- 
tion of the system along a complex trajectory - solution 
to the classical equations of motion analytically contin- 
ued to the complex values of coordinates and/or time 
0,01 ■ Given the complex trajectory, one calculates the 
tunneling probability at small values of the Planck con- 
stant h, 



(1) 



where F and A are the suppression exponent and pref- 
actor, respectively. 

Recently a new, intrinsically multidimensional, mech- 
anism of tunneling has been discovered 0, @] ■ It differs 
qualitatively from the well-studied case of direct tunnel- 
ing, where complex trajectories connect the in- and out- 
regions of the phase space. The new mechanism gener- 
ically occurs in the situation when the total energy of 
the system exceeds the height of the potential barrier 
separating the in- and out- states but, still, the process 
remains exponentially suppressed (dynamical tunneling). 
The complex trajectories in the new mechanism end up 
on a real unstable periodic orbit lying on the bound- 
ary between the in- and out- regions. We call this or- 
bit "sphaleron"; its instability implies that the tunnel- 
ing trajectories are also unstable. The above behavior 
of complex trajectories leads to the physical picture of 
tunneling as a two-step process 0, HI : formation of the 
sphaleron and its decay into the out-region. The latter 
step is not described by the tunneling trajectory; on the 
other hand, it does not involve exponential suppression, 
as the decay of the sphaleron proceeds classically. This 
tunneling mechanism is generic and has been found in 
several quantum mechanical [l^, 11 1 and field theoreti- 
cal {ill models. It is natural to call the new mechanism 
"sphaleron-driven" tunneling [l^ . 

Tunneling via unstable semiclassical solutions raises a 
number of issues. First, search for unstable trajectories 



is problematic from the numerical point of view. Second, 
even if one finds the tunneling trajectory which tends to 
the sphaleron as i ^ +oo, a problem remains to describe 
semiclassically the subsequent decay of the sphaleron or- 
bit. Yet another issue is the calculation of the prefactor 
A. In the case of direct tunneling this calculation involves 
the analysis of linear perturbations around the tunnel- 
ing trajectory. This procedure is not applicable in the 
sphaleron-driven case, when the perturbations destroy 
the tunneling solution completely. 

In this Letter we systematically develop a general 
method to solve the above problems. The idea is to in- 
troduce a constraint into the path integral for the tunnel- 
ing amplitude. This modifies the equations of motion in 
such a way that tunneling trajectories are pushed away 
from the sphaleron. The modification is governed by a 
regularization parameter e. At e > the semiclassical so- 
lutions interpolate between the in- and out-regions and 
are stable. They describe the whole two-stage process 
of sphaleron-driven tunneling. Integration over the con- 
straint corresponds to taking the limit e +0. The orig- 
inal unstable trajectory is recovered from the regularized 
solutions in this limit. We find expressions for the sup- 
pression exponent and prefactor of the sphaleron-driven 
tunneling probability in a closed form. Our method is 
based on the ideas of 

Our analysis reveals, among other things, one universal 
feature: the prefactor A in the probability of tunneling 
via unstable trajectories gets suppressed by an additional 
factor ft^/^ as compared to the case of direct tunneling. 
This, in particular, implies non-trivial properties of the 
transition between the direct and sphaleron-driven tun- 
neling regimes. 

Remarkably, the method of e-regularization can be 
used to deform real solutions corresponding to purely 
classical motion into trajectories describing tunneling 
0, [IBl ■ This makes the method efficient for finding and 
classifying complex solutions in the case of chaotic tun- 
neling JlO|. We stress, however, that the phenomenon of 
sphaleron-driven tunneling is unrelated to chaos and was 
observed both in chaotic and regular systems. 

While our approach is completely general, for concret- 
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ness, we illustrate it using a two-dimensional model with 
the Hamiltonian 

H^{pI+pI + Lu'y^)/2 + exp[-(x + y) V2] • (2) 

The model describes the motion of a particle in a po- 
tential valley extended along the cc-axis with quadratic 
confining potential in the y-direction. The valley is in- 
tersected at an angle by the potential barrier which in- 
troduces non-linear coupling between the degrees of free- 
dom. The process we are interested in is a penetration 
through the barrier of the particle which comes from the 
left in a fixed initial quantum state \E,Ey). The latter 
is characterized by the total energy E and the energy of 
oscillations in the y-direction, Ey = huj{n + 1/2), where 
n is the occupation number of y-oscillator. Note that we 
keep Ey fixed in the semiclassical limit h 0, so that n 
grows to infinity. It is shown in Q that, for given Ey, 
transmission through the barrier is a tunneling process 
for total energies E < Eb{Ey), while at > Eb{Ey) the 
transmission proceeds classically. The values of Eb are 
considerably higher than the height of the potential bar- 
rier, Eb{Ey) > Vmax — 1- The mechanism of transmis- 
sion changes from direct tunneling to sphaleron-driven 
tunneling when the total energy exceeds a certain criti- 
cal value Ec{Ey), where Vmax < Ec{Ey) < Eb{Ey). 

We start by reviewing the derivation of the formula 
([1]) in the ordinary regime of direct tunneling. One con- 
siders the amplitude of transition from the state \E, Ey) 
at the initial time moment t = ti to the state |x/) with 
definite coordinates beyond the barrier at the final mo- 
ment t — tf. Using the propagator in the coordinate 
representation, one writes. 



where 



A = J d^,(xf\e~'"^'/^\^,){K,\E,Ey 



(3) 



where At = tf — ti. The propagator is given by the 
semiclassical Van Vleck formula, 

1/2 

, (4) 



iS"^'(xi,xy)/S 



2'Kih 



det 



where S''^'(xi,xj) is the action evaluated on the classical 
solution going from at t = to at t = f/. It will 
be important for us that ([4]) can be derived [l3| from the 
path integral 



(x/|e-'^^*/''|x,) = 



x(t/)=x/ 

[dx(<)] e*^["(*)l/'^ . (5) 

We take the wave function of the initial state in the form 
{yii\E,Ey) — ilj{xi)'^{yi), where ^'(s^i) is the plane wave 
with the unit flux normalization, and '^(jji) is the semi- 
classical expression for the oscillator eigenf unction. Eval- 
uating the integral ([3]) in the saddle-point approximation 
one obtains, 



A 



2ttD^ 



^'i{S''' + B}/h+iTr/4: 



(6) 



B — XjXj, 



Di = Xiiiidet 



dy\/2Ey-iu^y'^ 




All the quantities in ^ are evaluated on the saddle tra- 
jectory satisfying the initial conditions 



= 2{E-Ey) 



yf+uj^y^^2Ey , 



(7) 



where x^ = x(ti). In deriving ^ we assumed that the 
saddle configuration is unique; this is indeed the case for 
the model ©. 

The amplitude (O is to be inserted into the formula 
for the tunneling probability. 



if , (8) 



lim 4- / dxflAf ^ [ dytx 



where in the second equality we cancelled the diver- 
gence of the integral originating from the infinite region 
— > +00 by writing J dxf = if - Note that this op- 
eration is legitimate only if the boundary condition 



as 



t 



f 



(9) 



is satisfied. One substitutes (O into ([S]) and evaluates 
the saddle-point integral over yf. One arrives at the ex- 
pression ID) with 

where D2 — 2 Im ^jj; . The saddle-point equation de- 
termines the remaining boundary condition for the tun- 
neling trajectory. 



Imyy = Imyj = . 



(11) 



Note that reality of the total energy together with ([5]), 
(fTTj) imply that i/ is also real. 

One solves numerically the classical equations of mo- 
tion with the boundary conditions JT]) , © , pT|) and dis- 
covers Q that tunneling solutions ^] satisfying all the 
conditions exist only at £' < Ec{Ey). Still, one can find 
solutions at Ec < E < Eb by imposing, instead of ([9]), 
the reality condition Imxf = Imxf = 0. But the re- 
sulting trajectories never come out from the interaction 
region: they end up oscillating on top of the potential 
barrier. These solutions att^ +00 become precisely the 
sphalerons we referred to above. The formula pO)) for the 
prefactor is not applicable in this case: the quantities Di , 
D2 entering into it describe linear response of the tun- 
neling solutions to small perturbations of the boundary 
conditions and are ill-defined due to the instability of the 
solutions. 
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To deal with this situation, we make the following 
steps. First, we introduce a functional Tint[x{t)] defined 
on classical paths. The choice for Ti„t[x(t)] is restricted 
by three requirements: (a) it must be real and positive- 
definite on real paths, (b) it must be finite on paths end- 
ing up in the out-region, (c) it must diverge on paths 
which stay forever in the interaction region. Overall, 
the functional Ti„t[x(t)] should roughly measure the time 
spent by the particle in the interaction region. For the 
model ^ the simplest choice is Ti„t[x(i)] ^ J dt /(x(t)), 
where the function /(x) > vanishes at x — > ±00. 

Second, we restrict the path integral (O to paths stay- 
ing fixed time in the interaction region, Ti„t[x(i)] = r. 
This eliminates unstable trajectories from the domain of 
integration. The full propagator is then recovered by in- 
tegrating over T. This program is realized by inserting 
the unity 



dT5{T,Mt)]-T) = jdrj ^e' 



(r-T„t[x])/n 



into ([5]) and changing the order of integration. We obtain. 



(x/|e 



-iHAt/h\ 



= / dT 



ide 



r/h 



+iQO ^'^'^ Jx(ti)=Xi 



,i(S[x]+ieT„.t[x])/?i 
(12) 



One observes that the integral over [dx] in (fT^ is the 
same as in ([5]) up to the substitution 



S i-^ Sc = S + ieTi, 



(13) 



Therefore, one can follow the steps leading from ([5| to 
H with the result 



A 



dT 



2Trh 



r/h 



^i(Sf+B,)/h+iTr/4 



2ttD 



Importantly, the semiclassical trajectory Xe(t) here is a 
solution to the equations of motion obtained from the 
modified action Se- By construction, it spends a finite 
time interval in the interaction region, and thus is stable. 
The integral over e is saturated by the saddle-point e(r) 
which is implicitly defined by the condition Tint[xg] = r. 
The latter follows from the relation d{Sf + Bf)/de = 
iTint\pi-(\- Note that the saddle-point value e(T) is not 
necessarily purely imaginary; one should be careful to 
pick the one with Ree(r) > in order to ensure the 
convergence of the path integral in ([T2|) . 

One proceeds by substituting the expressions for the 
amplitude and its complex conjugate into the tunneling 
probability and performing integration over the final co- 
ordinates. This leaves the integral over two interaction 
times, T and r', coming from the amplitude and the com- 
plex conjugate amplitude. It is convenient to change the 
integration variables to r_|_ = (r -f t')/2, t_ = r — t'. 



The integral over r_ is, again, saturated by the saddle 
point; one uses the formula d{Sf + — ieT)/dT = —ie 
in deriving the saddle-point condition. One obtains 



V 



UJX,J^~de/dT+ ^_2{Im{Sf+B,)-.r+)/n 



where the integral is performed over the real axis, and the 
function e(T_^) is defined by the folowing implicit relation 



r„t[Xe] -t- r„f [x_e] = 



(15) 



Note that the solution to (|15p is generically real for real 
T-f.: the equations of motion following from lead to 
x_e. = X*, which for real e implies that the l.h.s of (fTS]) 
is, indeed, real. 

So far, we did not refer to the particularities of the 
sphaleron-driven tunneling. The step where they be- 
come important is integration over r-|_. For both direct 
and sphaleron-driven tunneling, the integral (1141) is dom- 
inated by the point e(r-)_) = which corresponds to the 
original unregularized trajectory. In the standard case 
of direct tunneling the unregularized trajectory spends 
a finite time interval in the interaction region; one 
obtains the expressions (|10p for the suppression expo- 
nent and prefactor by the saddle-point integration over 
T+. In the sphaleron-driven case the interaction time 
corresponding to e = is infinite. Thus, the integral 
is saturated by the end-point of the integration inter- 
val, r+ +00. The tunneling probability in this case is 
determined by the behavior of the integrand in (fT4|) at 
T+ —f +00, that is, e —f +0. Performing the integration, 
one obtains the formula ((T|) with 



F = lim Fe 

£—+0 



A=h}^'^ lim ■ 



A, 



-47r 



d Rc Ti„t [x, 
dt 



, (16) 



where and are given by (|T0)) evaluated on the e- 
regularized solution x^. [A subtle point is that F^ = 
2{lm{Sf + B^) - eT+) = 2Im(S'[xe] + B^) is to be com- 
puted using the original action evaluated on the regu- 
larized solution.] Note an additional factor h}/"^ in the 
prefactor compared to the case of direct tunneling. 

Let us summarize our results. We have derived the 
following prescription for the calculation of the tunnel- 
ing probability in the sphaleron-driven case. First, one 
replaces the action of the system with the modified action 
St: P^ . where e > 0. Second, one finds the tunneling 
solution Xg of the modified equations of motion. This so- 
lution interpolates between the in- and out-regions and 
is stable. One evaluates its suppression exponent F^ and 
prefactor using the ordinary "direct tunneling" for- 
mulae. Third, one determines the suppression exponent 
and prefactor by taking the limit e -f according to 
P^ . Note that our prescription does not make use of 
any particular properties of the illustrative model ([2]); 
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it is applicable in a large class of models exhibiting the 
phenomenon of sphaleron-driven tunneling. 

We checked the method of e-regularization by applying 
it to the model ^ and comparing the semiclassical re- 
sults with the exact suppression exponent and prefactor 
extracted from the numerical solution to the Schrodinger 
equation. The latter quantities were obtained by fitting 
the dependence of the exact tunneling probability on h 
with the formula 

hlogVih)Ki-F + ^h\ogh + hlogA, (17) 

where A is independent of h. For completeness, we also 
performed the comparison in the regime of direct tunnel- 
ing. We set 7 = 1 for the direct tunneling and j — 2 
in the sphaleron-driven case. Figure [1] shows the de- 
pendences dl F{E) and A{E) at fixed Ey = 0.05. The 
semiclassical and exact quantum-mechanical calculations 
are in good agreement. 

Note that the quality of the fit ^T7} becomes worse as 
one approaches the transition point Ec between the two 
tunneling regimes. This is a manifestation of the break- 
down of the semiclassical approximation in the vicinity 
of this point. It is caused by the change in the depen- 
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FIG. 1: Comparison between the semiclassical (lines) and ex- 
act (points) results for the suppression exponent and prefactor 
in the cases of (a) direct and (b) sphaleron-driven tunneling. 
The comparison is performed in the model ([2]) with uj = 0.5, 
Ey — 0.05. The error bars represent the uncertainty of the fit 



dence of the prefactor on h. Indeed, A{E) diverges as 
E ^ Ec + 0, see Fig. [TJj, which contradicts the con- 
tinuity of the exact tunneling probability. However, the 
size of the vicinity where the semiclassical approximation 
breaks down vanishes in the limit S — > 0. 
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